arXiv: 1506.0821 Ivl [math.AP] 26Jun2015 


Exponential decay estimates for the stability of boundary 
layer solutions to Poisson-Nernst-Planck systems: one 

spatial dimension case 

Chia-Yu Hsieh * Tai-Chia Lin ^ 


Abstract 

With a small parameter e, Poisson-Nernst-Planck (PNP) systems over a finite one¬ 
dimensional (ID) spatial domain have steady state solutions, called ID boundary layer 
solutions, which profiles form boundary layers near boundary points and become flat 
in the interior domain as e approaches zero. For the stability of ID boundary layer 
solutions to (time-dependent) PNP systems, we estimate the solution of the perturbed 
problem with global electroneutrality. We prove that the norm of the solution of 
the perturbed problem decays exponentially (in time) with exponent independent of e if 
the coefficient of the Robin boundary condition of electrostatic potential has a suitable 
positive lower bound. The main difficulty is that the gradients of ID boundary layer 
solutions at boundary points may blow up as e tends to zero. The main idea of our 
argument is to transform the perturbed problem into another parabolic system with a 
new and useful energy law for the proof of the exponential decay estimate. 


1 Introduction 

The Poisson-Nernst-Planck (PNP) system, a well-known mathematical model for ion trans¬ 
port, plays a crucial role in the study of many physical and biological problems [D a El El El El 
[Hiiniiieiiniiii]. Such a model can be represented as 

Tlf ^ ' Jn •) Pt ^ ' Jp •) (^*^) 

Jn = -D^ , J, = -D, {vp + , ( 1 . 2 ) 

£A0 = — p + ZnCn — ZpCp , (1.3) 

for a: G > 0, where {n,p,(j)) depends on x and t, Q <Z is a bounded smooth domain 

N 

in R^,iV > 1, V = {dxi,-'' and A = the Laplacian. Physically, 0 is the 

i=i ' 
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electrostatic potential, n is the charge density of anions, p is the charge density of cations, 
p is the permanent (hxed) charge density in the domain, Zp are the valence of ions, e 
is the elementary charge, ks is the Boltzmann constant, T is temperature, Jn, Jp are the 

ionic flux densities and -D„, Dp are their diffusion coefficients. The parameter e related to the 

dielectric constant and the Debye length can be assumed as a small parameter tending to zero 
(cf. [21 [71 Us]). For simplicity, we only consider monovalent ions, that is, Zn = Zp = 1, and set 
e/ksT = 1, p = 0, Dn = Dp = 1. Besides, we rescale e and transform (ll.ip - (ll.3p into 

nt = -'V-Jn, pt = -V-Jp, (1.4) 

=-(Vn - nV0), Jp =-(Vp + pV0), (1.5) 

eA(j) = n — p, (1.6) 


for X G D, t > 0. 

Debye (diffuse) layers occur in ionic liquids near electrodes and have many applications 
in the helds of chemical physics and biophysics (cf. [ID]). To see Debye layers, solutions of 
(ll.4|) - (ll.6|l with boundary layers need to be investigated. For simplicity, the domain D is set as 
D = (—1,1) a one-dimensional interval in the whole paper. Then (II.4^ - 01.6j) can be denoted as 


{ nt = d^{nx - n(j)a:) , 

Pt = dcciPx + p4>x) , (1-7) 

£(j)xx = n-p, 

for X G (—l,l),t > 0. For the boundary conditions of fll.7p . we consider no-flux boundary 
conditions of n and p to describe the insulated domain boundaries, which are commonly used 
to study physical (biophysical) phenomena like the electric double layer and the ion trans¬ 
port through channels. Besides, we use Robin type boundary condition of cj) to represent the 
capacitance effect of physical systems (cf. [H] [13 [IB]) given as follows: 

{ nx - ncf)x = Px+ P<Px = 0, at x = ±1, 

0 + 7£0x = 0o( 1) at x = l, (1.8) 

0 - 7£0x = 0o(-l) at X = -1, 

where (/>o(l), 0o(~l) are constants, and > 0 is a constant depending on e. System fll.7p 
with fll.Sp has the conservation of total charges of the individual ions ndx = A, J^^pdx = 
B for t > 0, where A and B are positive constants (independent of t) representing total 
negative and positive charges, respectively. In most of the physical and biological systems, 
global electroneutrality holds true which means the total positive charge equal to the total 
negative charge. Consequently, we assume that A = B = mo > 0, i.e., global electroneutrality 
holds true in the whole paper. 

System fll.7p with fll.81) has a steady state solution (n,p, 0) = (n°,p°,'0) denoted as 

/= , (1-9) 

e'^dx e~'^dx 
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where tuq > 0 is a constant and "0 is the solution of the following equation, called charge- 
conserving Poisson-Boltzmann equation (cf. [IS]), with Robin type boundary conditions: 


^'4^ XX 

< 


mo mo e 

fli e~'^dx 


,'0(±1) ±7£'0x(± 1) = 0o(±l)- 


( 1 . 10 ) 


Note that — rP'ipx = pS + = 0 ^J^cl e'lpxx = — p^. Without loss of generality, we 

may assume 0o(l) = “0o(~l) > 0. From [12], we get the following results of boundary layer 
solutions of fll.lOp : 

Theorem A. (cf. [12]) Let 4j be the solution of fll.lOl) . Then 'ip is odd, i.e., ipix) = —'ijj{—x) for 
X G [—1,1], increasing in (—1,1), convex in (0,1), and concave in (—1, 0). Moreover, -0 satishes 

(i) Interior Estimate: 


10 (a:) I < 00 (1) 



^(l+x) 


+ e 


M 


(1-x) 


and lim 0 ix) 

e^O-l- 


0 for X G (—1,1) , 


where M = is a positive constant independent of e. 

(ii) Boundary Estimate: 

If 0 < lim = 7 < cx), then lim 0 (1) = 0* and lim v^0' (1) = -v/a — e 

£^ 0 + £^■ 0 + £^ 0 + \ 

where 0 < 0* < 0o (1) is uniquely determined by 0o (1) — 0* = 7\/« (e ^ — e ^ ). 


Theorem A implies that as e goes to zero, 0 have asymptotic behavior of boundary layer so 
we call (n,p, 0) = (n°,p°, 0) as a boundary layer solution of system (ll.7p with (II.8p . Note that 
(n°,p°) is represented in (11.91) . 

To get the stability of the boundary layer solution (n°,p°,0) to system fll.7l) with fll.Sp . we 
study the perturbed problem fll.lip with fll.l2p which comes from the assumption that system 
fll.7p with fll.Sp has solution 


(n,p,0) = (n0p°,0) -1- {n,p,Pj , 

where ^h, p, 0^ satishes the perturbed problem 

{ ^t TT'XX (^ 4^x')x {j^'4^x')x {j^(4x')xi 

Pt Pxx “1“ (P 0a;)a: d" (p0a;)a: "F ip4^x)x^ 
s^xx = n-p, for X G (-1,1) , t > 0 , 

with boundary conditions: 

jfix- n°03. - ni>x - n^x = Px + P°0x + Pi>x + P^x = 0 at x = ±1, 
I 0 ± 7£0a; = 0 at X = ±1 . 


( 1 . 11 ) 


( 1 . 12 ) 
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Here ft and p denote charge density perturbation of anions and cations, respectively. To let 
the global electroneutrality hold true, i.e., f_^ndx = A = B = J_^pdx > 0, we need to have 

fidx = J^^pdx = 0, which conies from condition (11.251) . Otherwise, ii A B and A, B > 0 
independent of £, then as £ approaches zero, steady state solution of fll.7p with fll.Sp becomes 
unbounded and far away from (n°,p°, '^) the boundary layer solution of fll.Tp with fll.Sp for the 
case of A = B = uiq > 0 (cf. [I2])- It seems impossible to get the stability of {rA,p^,'ilj) if 
condition (ll.25p fails and ndx = A ^ B = pdx (H, B > 0 independent of e) holds true. 
This motivates us to assume A = B = mo > 0 and (ll.25p in the whole paper. 

Conventionally, the stability of (ll.7p with e = 1 and the Dirichlet boundary condition for 0 
holds true because of lim^^oo ||h|| + IIpII ]^oo — 0 (cf. [3]) and the exponential decay estimate 

W^WlI + WpWlI < (cf. |1]) for f > 0, where the constant C and exponent A are positive. 

Here we study the stability of fll.7p with 0 < e <C 1 a small parameter tending to zero and 
the Robin boundary condition for cj) (see (II.8p ). It seems difficult to get the exponential decay 
estimate in LCnorm with exponent independent of e. The main difficulty is that the prohle of 
the solution ^|J has boundary layers near boundary points x = ±1, and blows up at boundary 
points a: = ±1 with order i.e., |'03;(±1)| ~ as e tends to zero (cf. [I2]). Instead 

of the LCnorm estimate, we prove the exponential decay estimate in i7“^-norm denoted as 
||h||^-i + ||p||^-i < for f > 0, where a is a positive constant independent of e, and Jq is 

a positive constant depending on the i7“Cnorm of initial data n\i^Q and p\^^Q■ 


1.1 Main Results 

To study system fll.lip with fll.l2p . we introduce the change of variables 

6 = n — p, f] = h + p, (1-13) 


where 6 is the gap between concentrations of positive (cations) and negative (anions) ions. 
Note that f^^Sdx = 0, i.e., f^^ndx = f^^pdx means the total positive charge equal to the 
total negative charge, and then global electroneutrality holds true. By (I1.13p . system (II.lip 
with (I1.12P becomes 

{ ^t ^XX (p 4^x')x (^P'4^x')x (^p4^x')xy 

Vt = Vxx - (^Vx)x - iSi>x)x - (50x)x, (1-14) 

^4^xx 


with boundary conditions: 

- pi)x - P^x = Px- d^^x - dipx - 50X = 0 at X = ±1, 
|0±7e0,j; = O at a: = ±1, 


where due to fll.131) . 


(5° = 


7° = , 


(1.15) 


(1.16) 


and {nP,p^) is dehned in fll.9p . By fll.9p . fll.lOl) and fll.161) . '0 satisfies 

eipxx = d^ for a:e(-l,l). (1.17) 
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Linear Stability 

To get linear stability, we consider the linearized problem of fll.ldp with fll.lSp as follows: 


^xx (h 4^x')x 

ht Vxx 4^x')x (^t/^a;)x) 


with boundary conditions: 


^x-V^^x - vi^x = Vx - ^°^x - = 0 at x = ±1, 


± 7£03; = 0 at x = ±1. 
By fll.181) and fll.lQp . it is obvious that 


d 

dt 


l_i dt 


fjdx = 0 for f > 0 , 


'-1 


and then 




6{x,t)dx = J 6o{x)dx, 
fj{x,t)dx = / fiQ{x)dx , 


'-1 


'-1 


for all f > 0, where 5o(2^) = 0) and fio{x) = fj^x, 0) are the initial data. Let 


and 


D{x,t) = J 5{s,t)ds 


H{x,t) = J r]{s,t)ds, 


(1,18) 


(1.19) 


( 1 . 20 ) 

( 1 , 21 ) 


( 1 , 22 ) 


(1,23) 


for X G (—1,1) and t > 0. Then by fll.20l) - fll.23p the boundary conditions of D and H at 
X = ±1 become the zero Dirichlet boundary condition 


D = H = 0 at X = ±1, for f > 0, 
if the initial data So{x) = S(x, 0) and 0) satisfy 


6o{x)dx = / fio{x)dx = 0 . 


(1,24) 


(1.25) 


'-1 


'-1 


The physical meaning of (ll.25p is: the total charge of solution {n,p) = (n° + h,p°+p) is same as 
that of (n°, p°) so the (global) electroneutrality n (x, t) dx = nP (x) dx = (x) dx = 

ill P t) holds true for f > 0. 
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We use fll.22p . fll.23p . and integrate equations of fll.lSp from —1 to x. Then we get 

Dt = (1.26) 

Ht = hfxx - , (1-27) 

which give the energy law of {D, H) expressed as follows: 


Ijtf + 


'-1 


'-1 


Dl + H‘^ + -r]^D^)dx-(l)^{-l,t) I iri^D + 5^H)dx (1.28) 




for f > 0 (see Theorem 12.11) . However, 03 ;(—l,f) = — j^^Ddx (see (I2.13p in Section |2]) 

and 7 £ ~ ^/e^ as e —)■ 0+ (see Theorem A) imply that 03;(—l,t) becomes extremely large as 

e approaches zero if the integral Ddx is away from zero. This makes fll.28p hard to be used 
for the estimates of D and H. 

To overcome such difficulty, we use the following transformation: 

1 


D{x,t) = D{x,t)— 


D{x, t)dx , 


^ J-i 
1 

H{x,t) = H{x,t) — - H{x,t)dx, 

2 J-i 

which can be denoted as D{x,t) = D{x,t) — d(t) and H{x,t) = H{x,t) — h(t), where 

-1 


h{t) = 


d(t) = 2 / ’ 


'-1 


H{x, t)dx. 




Then we prove the following result for fll.26p and fll.27p . 

Theorem 1.1. Suppose l[1.25\} . 0 < < 'fmax < oo, and ^ > (d+^ 7 maH j-3)xo0o(i) true, 

where M = J > Q, Kq = sup ^ ^ Imax > Q is a constant independent of 

" 0<|y|<<Ao(l) 

£. Then 


fL 

dt 


(T)2 + H^)dx + d^ + h^ 


'-1 



C / 

At 

U 


D'i + \hI) - lm„. 




(1.29) 




1 — 77000 ( 1 )^ > ond e depends only on mo 


for t > Q and Q < e < e, where mo,e = mo 
and 0 o(l)- 

Remark 1.1. Constant Kq is defined by Kq = sup ^ ^ which depends only on 

o<hl<0o(i) 

'^e 

00 (1) approaches to one as 0o (1) tends to zero. Hence for any 7 > 0, z/ lim —A = the 

e^0+ y'e 

hypothesis ^ > (d+^T'maH j-3)i<'o</)o(i) fulfilled z/ 0 o(l) > 0 zs sufficiently small or mo > 0 

is sufficiently large for all small e, where M = > 0. 
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From Theorem 11.11 we obtain the following estimates. 


Corollary 1.2. // 0 < 7 ^ < '^raax < 00 and ^ 
positive constant a > 0, independent of e, such that 


then there exists a 


1 

2 


+ H^)dx + d^ + h‘^ < Joe 


—at 


f-1 


(1.30) 


for t > 0, where 


h = 


(T)2 + H^)dx + d^ + h^ 




t=o 


(1,31) 


Because of D = D + d and H = H + h, Corollary 11.21 gives ||hl||^2 + ||hf||L2 < IqC i.e., 
||h||^-i + ||p||^-i < for f > 0, which implies the linear stability of fll.7l) with fll.81) in 

77“^-norm, where Jq is a positive constant depending only on the initial data. Here we have 
used the equivalence between ||ll||j ;^2 + ||hf|lL 2 and ||h||^-i + ||p||j|^-i (see Appendix I). 

To get the TT^^-norm estimate, we hrstly transform the linear part of the perturbed prob¬ 
lem fll.lll) with fll.121) (i.e., fll.lSp . fll.22p and fll.23p i into a coupled system of linear parabolic 
equations of {D, H) denoted as fll.26p and fll.27p with zero Dirichlet boundary condition fll.241) . 
To preserve the global electroneutrality, we assume that the total charge density perturbation 
is zero for anions and cations, i.e., the initial data satishes (11.251) . which implies boundary con¬ 
dition fll.241) . Then we hnd the associated energy law fll.281) (proved in Theorem 12.11) but the 
coefficient of the last integral of fll.281) still blow up as e tends to zero if the integral D dx 
is away from zero (see fl2.13p '). This motivates us to decompose {D, H) into (Z), ZZ) and (d, h), 
where D = D — d, H = H — h, d = \ D dx and h = Z H dx. Then we derive fl2.19p (see 
Theorem 12.21) as the energy law of [D,H,d,h) to prove Theorem 11.11 and Corollary 11.21 which 
imply ||Z 1||^2 + ||ZZ ||^2 < Zoe“"* and hence ||h||^-i -|- ||p||^-i < Zoe“"* for f > 0 , which gives 
the linear stability of fll.Zp with fll.8p in ZZ“^-norm under global electroneutrality. Here a is 
a positive constant independent of e and Iq is a positive constant depending on the L^-norm 
of initial data and i.e., the ZZ“^-norm of initial data and Note that 

||Zi)||j ;^2 + ||ZZ||j ;^2 is equivalent to ||h||j|^-i -|- ||p||j:^-i (see Appendix I) and ye > 0 is assumed to 
have a suitable positive lower bound (see Theorem ll.ip . which makes the last three terms of 
fl2.40p together become nonpositive so fll.29p holds true. Such an assumption of 7 ^ is also used 
to study nonlinear system fll.141) with boundary condition fll.l5p . 

Nonlinear Stability 

For nonlinear stability, we may generalize the idea of linear stability to study {6, fj, (f) the 
solution of nonlinear system fll.l4p with boundary condition fll.lbp . The main difficulty is to 
control the nonlinear terms rjcfx and of system fll.l4p . Here we assume that the initial data 
satishes fll.33p . which implies that the right side of fll.32p becomes negative (see Theorem II.3p . 
Consequently, fll.32p is useful to show fll.34p (see Corollary II.4p and get the nonlinear stability 
of (?7,°,p°,'0) to system fll.Zp with fll.8p . 

Now we state results for nonlinear stability as follows: 
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Theorem 1.3. Under the same hypotheses as in Theorem \1.1[ suppose furthermore that the 
initial data {n,p)\^^Q = (rioPo) (~1) 1) satisfies no,po G (—1,1) and 

no{x) = n^{x) + h{x, 0) > 0 , Po{x) = p^{x) + p{x, 0) > 0 for x G (—1,1). 


Then 


d 

dt 


1 

2 



{D^ + H^)dx + d^ + h^ 



moE 


Hi I dx — —m'n 


Ae 


0 ,£ 


'-1 


D^dx 

(1.32) 


for t > 0 and 0 < e < e', where = mo 1 — 2i^o0o(l)^ > depends only on mo 

and </>o(l)- Moreover, if the initial data satisfies 


lo := 


+ H‘^)dx + d? + h^ 




< 


OmoS 


(1.33) 


J t=o 


for some 0 < 6* < 1, then ^ —— > |(1 — 6^) > 0 for all t >0. 


moe 


Besides, from Theorem 11.31 we get 

Corollary 1.4. Under the same hypotheses of Theorem 11.31 if U.33\) holds true, then there 
exists a positive constant a' > 0, independent of e, such that 


1 


( Z }2 ^ ^ ^2 < 


-a’t 


(1,34) 


'-1 


and 


l|0|li| + l|J^llii< V-“'‘ (1.35) 

for t > 0 and 0 < £ < U. 

Due to D = D + d and H = H + h, fll.34p may imply fll.351) and show that the upper 
bound of ||D ||^2 + ||-f^|li 2 being equivalent to ||h||j^-i + ||p||j|^-i (see Appendix I) approaches 
zero exponentially with exponent independent of e as t goes to infinity. This represents the 
exponential decay estimate (to e) of ||h||j|^-i + ||p||j(/-i and gives the nonlinear stability of fll.7p 
with fll.81) in norm. 

For nonlinear stability, we use the same idea of linear stability to study {S, fj, fi) the solution 
of nonlinear system fll.lip with boundary condition fll.151) . The main difficulty is to control 
the extra nonlinear terms fjcfx and dfix of system fll.lip . Here we assume that the initial data 
(^,p)lt=o = (^o,Po) in (-1,1) satishes no,po e (-1,1), no{x) = ?7,°(x) +h(x, 0) > 0, po{x) = 
p^{x) +p{x, 0) > 0 for x G (—1,1) and fll.331) which expresses the smallness of ||h||j:^-i + ||p||j|^-i 
at t = 0. Then we use fl3.4l) the energy law of [D, H, d, h) to show ||D ||^2 + ||hf||i 2 < for 

t > 0 (see Theorem 11.31 and Corollary 11.41) . i.e., ||h||^-i + ||p||j:^-i < Ioe~°‘'^ for t > 0, where a' 
is a positive constant independent of £, and constant /q > 0 comes from the L^-norm of initial 
data D\^^q and i.e., the of initial data n\^^o and satisfying fll.331) . Note 

that condition no{x),po{x) > 0 for x G (—1,1) implies n{x,t),p{x,t) > 0 (see Proposition 13.2p 








































and ri(x,t) = n{x,t) + p{x,t) > 0 for x G (—> 0, which implies J^^rjD^dx>0 a 
crucial inequality for the use of (13.dh to prove Theorem 11.31 In physical point of view, the 
nonnegativeness of no and po is reasonable because no and po are concentrations of anions and 
cations, respectively, at the initial time t = 0. 

The rest of this paper is organized as follows: For linear stability, we prove Theorem 11.11 
and Corollary 11.21 in Section [H In Section [3l the proofs of Theorem 11.31 and Corollary 11.41 are 
provided for nonlinear stability. 

2 Proof of linear stability 

In this section, we study (ll.lSjl with (11.1911 . which is the linearized problem of (11.1411 with 
(11.1511 . We derive the energy law of {D, H) as follows: 

Theorem 2.1. Let {S,fi,(l)) be the solution of U.18\) with boundary conditions U.19\) . If U.25\) 
holds true, then {D, H) satisfies 

\jtj + = (^Dl + H^ + ^p^D^yx-^,{-l,t) j {p^D + 5^H)dx (2.1) 

for t > 0. 

Proof. For equation (ll.26p . we multiply it by H, integrate it from —1 to 1, and use integration 
by parts. Then 


1 d 

2 dt 






dx. 


( 2 . 2 ) 


Here we have used the fact that H = 0 at a: = ±1 from (11.2411 . On the other hand, we integrate 
£(l>xx = <5 the Poisson equation of (ll.lSp from —1 to x. Then 


e{^,,{x,t) - ^,,{-l,t)) = e j ^yy{y,t)dy 


= / hy,t)dy 


= D{x,t), 


which gives 

0„(a:, t) = ^D{x, t) + 0„(-l, t) 
for X G (—1,1), t > 0. Consequently, 


p^cfxDdx = - I p^D'^dx + (j)x{,—I,t) I p^Ddx. 


'-1 


'-1 


'-1 


(2.3) 


(2.4) 
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For equation fll.27p . we multiply it by H, integrate it from —1 to 1, and use integration by 
parts. Then 


/ H^dx = 
2dt 






(2,5) 


Here we have used the fact that if = 0 at x = ±1 from fll.24p . Moreover, we use fll.l7p . fl2.3p 
and integration by parts to get 


((5°0,ii + xlj.D,H)dx = 


'-1 


'-1 


“1“ ) Hdx 


= J e{'ifj:^(j)x)xHdx 

= - j eil)^^a^H^dx 

= -j ^|JxDH^dx - [ i!xHxdx 

/•i 


'-1 


i/jxDHxdx + (j)x{—l,t) / d^Hdx. 


( 2 , 6 ) 


'-1 


'-1 


Therefore, the proof of fl2.ip . i.e.. Theorem 12.II is complete by combining fl2.2p . fl2.4p . fl2.5p and 

fiTHD . □ 

In order to use Theorem 12.11 for the proof of the linear stability of {n^, , ip), we need to 

consider 03;(—l,f) the gradient estimate at the boundary point x = —1. Notice that cp satishes 


eipxx = d for X e (-1,1), 


with Robin boundary condition 


' ± = 0 at X = ±1, 


(2,7) 


( 2 , 8 ) 


for each t > 0. Fix i > 0 arbitrarily. Then we integrate both sides of equation fl2.7p in x over 
the interval (—1,1), and get 


e((px(l,t) - <Px(-l,i)) = £ / ^xxdx= / Wx = ii(l,t) = 0 , 


'-1 


'-1 


which implies 


Here we have used fl 1.221) and fll.24p . Thus fl2.8p and 02.91) give 

= -7£0x(l,i) = -7£0x(-l,t) = 


(2.9) 


( 2 . 10 ) 
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By dMl) and fICTD . we have 


20 ^(-l,t)= / {x^^)^dx 


'-1 


= / {(t)x + X(t)xx)dx 


'-1 


= - 20 (-l,t)+ / x^xxdx 




= -27e'#’r(-li<) + / 


( 2 , 11 ) 


'-1 


Furthermore, we use fll. 22 p (which gives = S), fl2.7p and integration by parts to get 


»i 


X(j)xxdx = - 


'-1 


£ J-i 

-f 

£ J-i 

1 /•' 


xSdx 
xD^dx 
Ddx. 


'-1 


Consequently, ( 12 . 111 ) and fl 2 . 12 p imply 

= 


Ddx. 


( 2 . 12 ) 


(2.13) 


2(1 + 7e)^ J-I 

Note that ~ 7/^7 as £ ^ 0+ (see Theorem A). Thus as £ approaches zero, 02 ;(—l,t) 

becomes extremely large if the integral Ddx is away from zero. This makes (12.ip hard to 
be used for the proof of the linear stability of {'nP, , ip). 

To overcome the difficulty, we transform D and H into D and H by truncating the average 
of D and H, respectively: 

1 D 

Ddx, 


D = D-- 


H = H - 




Hdx, 


1-1 


and 


d = - 


h = 


Ddx , 


-1 

rl 


Hdx , 


(2.14) 

(2.15) 

(2.16) 
(2.17) 


'-1 


where d and h are the average of D and H at time t, respectively. Note that D = D — d and 
H = H — h satisfy 


Ddx = / Hdx = 0 , 


(2,18) 


'-1 


'-1 


11 


















Dx = Dx = S and Hx; = Hx = r; for x G (—1,1) and t > 0. 

Now we state the energy law for (Z), ZZ, d, h) as follows: 

Theorem 2.2. Under the same assumption as Theorem \2.1[ we have 


for t > 0. 

Proof. Note that 


d 1 
dt 2 



+ H^)dx + d^ + h^ 


'-1 


Di + Hi + 


dx — 


2mo7gd^ 

(1 + le)£ 


(1 + 27 g)d 

(1 + le)£ 


rj^Ddx + 





6^Hdx 


'-1 


D‘^dx= / {D-dfdx= / D‘^dx-2d^, 


'-1 


'-1 


'-1 


j H^dx = j {H-hfdx = j H^dx- 


2H 


(2.19) 


( 2 . 20 ) 

( 2 . 21 ) 


By fl2.20p and fl2.2ip . equation fl2.ip becomes 


d 

dt 


+ H^)dx + d^ + h^ 


'-1 




( 2 . 22 ) 



Dl + Hi + -r]°D^)dx-(l)x{-l,t) I {v''D + 5^H)dx. 






Then we put fl2.13p into fl2.22p and get 


A 

dt 


(Zl2 + H^^dx + d^ + h^ 


'-1 


^ - 1 - 

Dl + Hi + -ri%D + dY)dx + 



d 


'-1 

'Dl + Hi + -r]^D^)dx- ^ 


(1 + 7e)^ J-1 


{7f{D + d) + (5°(ZZ + h))dx 





£ 


rf{2dD + d^)dx 


'-1 


d 


(1 + 7e)^ J-1 


+ 5^H)dx + 


d2 



'-1 


Di + Hi + -r]^D^)dx 

£ Z (l + 7.)e 


(1 + le)£ J-1 

2mo7£d^ 


rj^dx + 


dh 


(1 + 7£)e J-i 


6°dx 


{l + 2%)d 

(1 + 1 e)£ J -1 


rj^Ddx + 


d 


6^Hdx, 


>-i 
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which gives fl2.19p and complete the proof of Theorem 12.21 Here the last equality uses the fact 
that d and h are independent of x, and 



2mo, 



0 . 


□ 


2.1 Proof of Theorem 11.11 


In order to use Theorem 12.21 for the proof of stability, we need to estimate the integral terms 


of (1TT9D involving r/O and 5°, where r]° = 5° = ^ 




solution of fll.lOp satisfying 


/-i- 






and 'll; is the 


S'ipxx = 5^ for X e (-1,1). 


(2.23) 


By Theorem A, 


(a:)| < (/>o(l) e 




(2.24) 


for all X G [—1,1] and e > 0, where M = independent of e and 7 . Then we claim that 

< 2 + KoM^)^, (2.25) 

where Kq is a positive constant defined as follows: 

|e^ - 1 | 

Kq = sup —j—— > 0 . 

o<hi<<^o(i) \y\ 


Note that Kq only depends on 0o(l). By Theorem A, -0 is odd, increasing, and |'^(x)| < 0o(l) 
for X G [—1,1], which implies 



and 

p(x) < / 1 X > 0, 

- \ 1, X < 0. 

And then we may use fl2.24l) to get fl2.25l) . Now we claim that 

|?7° - mo - |7° - mol I < moA'o0o(l)^. 


We divide the domain interval [—1,1] into two parts as follows: 


A = {x G [—1,1] : 7 ° > mo} 


(2.26) 


(2.27) 


(2.28) 
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and 


-B = {x e [—1,1] : ?7° < mo} . 


Then we get 


rf — mo — \'rf — mo| I = 0 on A, 


and 


\ri^ — mo — \ri^ — mo| | = 2 (mo — rj^) on B. 
On B, by the dehnition of r]^ and fl2.26p . we get 


mo — ff = mo 

< mo 

< mo 


1 - ^-(e^ + e-^) 


1 - 


1 - 




< mo-ft'o0o(l) 


1 + i^o0o(l)^ 


2M' 


Here we have used the fact that + e ^ >2 and (12.251) . Therefore, we complete the proof of 
(12.281) . Moreover, by (I2.28p and the fact that = 2mo, we have 


|?7° - mo| dx < 2mo/i'o0o(l)^- 


(2.29) 




To get the gradient estimate of ip, we multiply fl2.23p by ipx and integrate it over (—l,a;). 
Then 


= TT^ (e’^ + e +0^, for x G (-1,1) , 


(2.30) 


where is a constant depending on e. Taking the value at x = 0 for both sides of fl2.30p . we 
have 


e ,2/n^ 2mo 


c, = -Mq) - 

^ 1:^1 

Integrate both sides of fl2.30p and by fl2.31l) . we get 


(2.31) 


£ 

2 


= 2mo 1 


'-1 


f-1' 


+ #x( 0 ). 
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Then using fl2.25p . we have 


e 

2 


< 2 mo 1 


'-1 


1 + i^o 0 o(l)^ 


+ #x( 0 ) 


< mo/^'o 0 o(l)^ + #x( 0 )- 


(2.32) 


By Theorem A, is increasing on [—1,1], convex in (0,1), and concave in (—1, 0), which implies 
i^x > l/^x(0) > 0. By the mean value theorem and (I2.24p . 


' 0 x(O) < 


i’ ( 1 ) “ !*(•)) 


= 2^12 


<20o(l)(e + e 2 v? 


Therefore, fl2.32p implies 
/•i 


In In/ 3M M \ 2 

I ^l^l< 2moKo(j)o{l)j^ + 40o(l) (^e + e 


And we have 


(2.33) 


(2.34) 


J ^rnoKo(poil)j^ 

for 0 < e < where Eq is a positive constant depending only on uiq and 0 o(l)- 
Moreover, we use fl2.23p and integration by parts to get 


(2.35) 


d^Hdx = £ ^l>xxHdx = —£ ipxHxdx. 




'-1 




The boundary integral of the last equality is zero because 'ipx is even and H{—1) = H{1) = —h. 
Hence by fl2.35p and Holder’s inequality, we have 


, ^ , f S^Hdx 

< 

(1 + 7 .)^ 7-1 



d 


< 


< 


l + 7e J- 

d 


1 + 7£ \J-1 

1 rl 


xfjxHxdx 

1 

1 \ 

^Idx] / Hldx] 


< - 


j-i 

, H^dx + , , , 

/_! " 2(1 + 7,)7_i 


d^ 




i^ldx 


1 - 3d^ 1 


(2.36) 


Here we have used the fact that 7 ^ > 0. 
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Now we want to estimate the integral 


Dne to Ddx = 0, 




as follows: 


rf^Ddx = — mQ)Ddx , 


1-1 


'-1 


so we may use the Holder’s inequality and ab < + ^6^, Va, 6 > 0, to get 

(1 + 2'y^)d 


(1 + J-1 

(1 + 2'yi;)d 


rfDdx 


{rf — mo)Ddx 


'-1 


< 


(1 + 2'y^)d 


(1 + 7^)^ 


1 Z' - \ 

\ri^ — mo\dx] ij {r]^ — mo\D^dx] 


<- \r]^ — mo\D^dx + 


(1 + 


£J-i' 4 ( 1 + 7 £)% 7 _i 

For the integral \v^ — mo\dx, we use (I2.29p 


| 7 ° — mo|dx. 




| 7 ° -mo\dx< 2 moi^o 0 o(l) 


M’ 


to get 


(1 + 


4(1+7j)% y_, 
Hence fl2.37p becomes 


\ri^ — mo I dx < 


(1 + 2%)d 


(1 + 75 )^ 7-1 


rf^Ddx 


^ 1 


£ J-i 


\'rf — mo\D^dx + 


(1 + 2%)^d‘^ moiFo0o(l) 

2(1 + 7 ,) 2 v ^ M • 

(1 + 2'^efd? moiFo0o(l) 


Furthermore, we may use fl2.28l) 


|7° - mo - |7° -"^oll < moi^o0o(l)^ 


2(1+7,)2v^ M 




and (12.381) to get 


£ J-i 


(1 + 27e)'^ f 


yfDdx < — J rf^D^ 


(1 + 7e)^ J-1 


(1 + 2%)d 


(2.37) 


(2.38) 


(1 + de)^ J-1 


rf*Ddx 


< -- 


(1 + 2%)^d‘^ mo/<'o 0 o(l) 


/_7" ^2(1 +7.)^7/5 M 


1 /" r I t 0 I 0 hi n 2 I (1 + 27 e)^(i^ moi^o<(’o(l) 

-JJmo+ir, -™„-|7 -m,\)]D + M 


< —mo 

e 


1 - i^o 0 o(l) 


VI 

M 


»i 


(1 + 2^^)^^^ mo/Fo 0 o(l) 


/ h 2 _ 

1-1 2(1 + 'Ze)V^ 


M 


(2.39) 
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Again, here we have used 7 ^ > 0. And we can choose £1 > 0 depending only on mo and 0o(l) 
such that 


for 0 < £ < £ 1 - We substitute fl2.36p and fl2.39p into fl2.19p . Then 


A 

dt 


{D^ + H^)dx + d^ + h^ 


'-1 


< - 


J-1 

1 

--mo 


Di + -Hi 


(2.40) 


1 - A'„,#.„(i) A] /‘b 2 


'-1 


(1 + 2'^^fd‘^ moA'o0o(l) 


2(1 + le)\/£ 

2mo7£;d^ 


+ 


M 

3d^ 


(1 + 7 £)£ 2 ( 1 + 7 ^) 


moA'o^o (1) 




Recall that 
(1136]); 

(EH: 


d 


(1 + 7e)^ J-1 


6^Hdx 


< 


Hzdx + 


3^2 


'-1 


2(1+7.) 


moA'o^o ( 1 ) 


My/e' 


dt 


i It (fl" + H^-)dx + <i= + h^\ = - /y (Dl + HI + \i>D^-)dx - 

- /t d'Xdx + /y 


2mo'ye<fi 


In order to get fll.29p . we need the nonpositiveness of the last three terms of fl2.40p together, 

which may hold true by assuming J^^)jCo(/)o(i) ^ ^g pg^yg used the fact that 

7 . + Imax- Therefore, we set e = minjeoAi} (^oAi > 0 only depends on mo and 0o(l)) and 
complete the proof of Theorem 11.11 


2.2 Proof of Corollary 11.2 


Recall that since D{±l,t) = H{±l,t) = 0 and JJ^Ddx = ji^Hdx = 0, we may use the 
Poincare’s inequality to get 

II-D||l 2 < c||zi2;|ii,2, ||-D||l 2 < cii/ixIUi, 

ll^llLE<c||hf.|U2, <C||^.IU2, 

where C is a positive constant from the Poincare’s inequality. Hence by the Holder inequality, 
we have 


\d\ = 


D 


'-1 


< -\\D\ 


Ll 


rl \ 1/2 

b 


^/ 2 , 


T>l 


L 2 ) 
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which implies |(i| < since Da; = D^- Similarly, \h\ < • Hence fll.29p 

implies 


d 

dt 


{D^ + H^)dx + d^ + h^ 


'-1 


< —a 


{D^ + H^)dx + d^ + h‘^ 


1-1 


for some a > 0 depending only on the constants from the Poincare’s inequality. Therefore, we 
obtain fll.301) and complete the proof. 


3 Proof of nonlinear stability 


To get nonlinear stability, we generalize the idea of linear stability to study (5, f], 0) the 
solution of nonlinear system fll.ldp with boundary condition fll.lSp . which also has conservation 
laws as follows: 

■qdx = 0 for t > 0 . (3.1) 


d r Ux=- 

dt dt j_i 


(3.2) 


As for linear stability, we assume the initial data of {6, fj) satisfying 

j 5{x,0)dx = J fj{x,0)dx = 0, 

which is same as fll.25p . Moreover, as for Theorem 11.11 we set 

/ X nx 

6{s,t)dS; H{x;t) = J fi{s,t)ds for a; G (—1,1), f > 0 , 

D = D — d and H = H — h, where d = | Ddx and h = ^ Hdx. 

To control the nonlinear terms 7)03, and dcpa of system fll.ldp . we assume that the initial 
data satisfies n {x, 0) = nP (x) + h (x, 0), p (x, 0) = (x) + p (x, 0) > 0 for x G (—1,1) and 

fll.33p . which implies that the right side of fll.32p becomes negative (see Theorem II. 3p . To 
prove Theorem II.3[ we first derive energy laws as for linear stability in Section [2l Such energy 
laws are represented as follows: 

Theorem 3.1. If {S, fj, 0) is a solution of with boundary condition U.15\) . and the initial 

data satisfy 113. Then we have 

^ ^ ^ {D^ + H^)dx 


2 dt 


(3.3) 


'-1 


Di + Hf + —r]^D^ + —r]D^)dx - M-ht) / {v^D + 6^H)dx, 
f-I ^ 2s 2s / J_i 


and 


fL 

dt 


{D^ + H‘^)dx + d^ + h^ 


'-i 


-1 


2mo%d‘^ 


Dl + Hl + ^p^D^ + ^pD^dx , 

^ 2e 2e ) (1 + 75 )^ 

“1 a /■! 


(3.4) 


(1 + 2'^;f)d 


d 


d 


, , , n^Ddx — - j Ha-Ddx + , , , 

(1+7.)£7-i ej-i {l + %)sj_. 


d^Hdx, 
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where rj = n + p, {n,p,(j)) is the corresponding solution of the PNP system ( [i. 7| j-( l77^) with 
initial data satisfying n (x, 0) = (x) + n (x, 0), p (x, 0) = p^ (x) + p (x, 0) > 0 for x G (—1,1). 

Proof. The proof of Theorem 13.11 is similar to those of Theorem 12.11 and Theorem 12.21 the 
difference is to deal with the nonlinear terms. By integrating the equations for 6 and rj in 
(11.141) from —1 to X, we obtain 


Dt Dxx h f^xHx (fxHx.) 

Hf Hxx h (fx f^X^X (fxldx- 


Multiply 03.51) by D, integrate it from —1 to 1, and do integration by parts, then 


1 d 

2 dt 





+ ) dx. 


(3.5) 

(3.6) 


(3.7) 


Here we have used the fact that D{—1) = D{1) = 0. On the other hand, we multiply 03.61) by 
H, integrate it from —1 to 1, and do integration by parts, then 

~ j H^dx = - j + 6^^xH + fjxDxH + ^xDxHyx. (3.8) 


Similarly, we have used the fact that H{—1) = H{1) = 0. By the same argument in Theorem 
EH we have 


p^cfxD + f^xHxD + d^^xH + fjxDxHjdx 


^ J-i 

For the nonlinear terms, we have 


dx + (j)x{—^,t) J {rfD + 5^H)dx. 


(3.9) 


'-1 


(fx^xD T f^xhdtxdPj dx J (fxi^dDPl^xdx 


rDHdx 


'-1 


1 


1 

£ 


6DHdx 


-1 

rl 


DxDHdx 




2e 


D^Hxdx 


'-1 

/•i 


= — / fjD^dx. 


2s 


'-1 


(3.10) 
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Notice that r] = + fj. Combining fl3.7p - fl3.10p . we obtain fl3.3p . Now we use the relations 

between D, H, D, H, d and h, and fl3.3p . Then 


d 

dt 


p2 ^ + d^ + h^ 


1-1 




Ol + H^ + -^ri°D^ + ^riD^)dx-M-i,t) / (rfD + S'‘H)dx 


(3.11) 


-1 

pi + Hl + + dY + 

-0,(-l,t) j {rj^D + d) + 6YH + h))dx 


'-1 


As for fl2.13p . we have 




1 

2(1 + 7e)e 



d 

(1 + 7s)£ ’ 


(3.12) 


because the Poisson’s equation of 0 here is the same as that of the linearized problem. Then 
fl3.11l) becomes 


d 

dt 


{D^ + H‘^)dx + d^ + h^ 


'-1 


Dl + Hl + —rp{D + dY ^viD + dY) dx 
1 

O/n , , rO/ 


d 


(1 -0 7 e )£ 


{r]YD + d) + 6YH + h))dx 


-1 


Dl + H! + ^^rfD^ + ^riD'^]dx 


u 


^ / v'^Ddx - - 

(1 + 72)^ J-i ^ 


2 mo7gd^ 
(1 + %)£ 
d 


d - - 

- I H^Ddx + — -— 

-1 (1 + le)£ 


S^Hdx. 


'-1 


(3.13) 


Here we have used rj = ■ Then we complete the proof of Theorem 13.11 □ 

Proof of Theorem 11.31 

In order to use fl3.4p for the proof of Theorem 11.31 we need to estimate terms in the right- 
hand side of fl3.4p . For the term jl^rjD^dx > 0, we need the nonnegative sign of r; = u -|- p 
which may come from the following result: 

Proposition 3.2. Let {n,p, 0) be the solution of /jl. with initial data no,Po ^ (“1; 1) 

and UqjPo > 0 for x G (—1,1). Then n,p > 0 for x G (—1,1), t > 0. 

The proof of Proposition 13.21 is standard and is given in Appendix II. A similar proof can be 
found in [3]. Proposition 13.21 implies p{x,t) = n{x,t) + p{x,t) > 0 for a; G (—l,l),f > 0. 
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Here we assume initial data no(x) = n (x, 0) = vP {x) + n (x, 0) > 0 and Po{x) = p{x,0) = 
(x) + p(x,0) >0 for X e (— 1 , 1 ). 

Now we need to deal with the last three integrals in the right-hand side of 03.41) . For the 
integral J^-^^p^Ddx, similar to 02.37p . we have 


(l + 27e)d /•' 

(1 + 7e)^ 7-1 


p^Ddx 


< — / \p^ - mo\ D'^dx + 


(1 + 2p,fd^ 




\p^ — mo\dx. (3.14) 


2(1 + 7s)% 7-1 

As for 02.37h . we use ab < Va, & > 0, instead of ab < Va, 6 > 0. We also use 

02.291) to estimate “ ''^ol dx and the fact 7 ^ > 0 to get 


(1 -|- 2'yP)d 


(1 Pe)£ 7-1 


p^Ddx 


< 


< 


2e 

1 


'-1 


|?7° — mo I D^dx + 


2e 


moJFo0o(l)(l 2p^Yd‘^ 

(1 + %yMy/e 

jri m\Udx+ + ■ 


(3.15) 


Furthermore, we use 02.281) to get 


—^ / p^D'^dx 


< -- 


< -- 


< 


2e 

27 

7, 

mp 

' 2e 


'-1 


p^D'^dx + 


'-1 


(1 + 2ps)d /•' 

(1 + 'le)^ 7-1 

(1 -|- 2%)d 


p^Ddx 


(1 + 'Ps)^ 7-1 


p^Ddx 


f ^ 0 I 0 hi f^2j I ^0-^000(1)(1 + 27£)^(i^ 

^[mo+(7 -mo-l7 -mp\)\Ddx+ - ^ - 


1 - Kp(t)p{l) — 


1 '^ ^ 2 ^^ ^ moA:o0o(l)(l + 27£)^d^ 


By 02.36p . we have 


d 


(1 + 76:)^ 7-1 


6^Hdx 




H^dx + 


(1 + Pe)M^ 

3d^ 


mpKp(j)o (1) 


(3.16) 


(3.17) 


r-i 2(1-|-7^) ^ ^ My/e 

for 0 <7 < 7o) where £0 comes from 02.35P and depends only on mp and 0o(l)- Holder’s and 
Young’s inequalities give 


d 


H^Ddx 


'-1 


< 


d 


Hldx 


mp 


-1 

1 



D^dx 


-1 


< / D dx + 


d^ 


4e 


’-1 


mpE J_i 


Hidx. 


(3.18) 


Recall 03.4p : 


A 

dt 


(F)2 + H‘^)dx + d/ + h^ 


'-1 



'-1 


" ^ ^ 2e^ 2e^ ) (1 + p/e 


/ + ‘2le)d r 

(1 +%)£ 7-1 


d r - - 

p^Ddx - J H^Ddx + 


d 


(1 + le)^ 7-1 


6^Hdx. 
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Substituting fl3.16p - fl3.18p into fl3.4p . we get 


d 

dt 

< - 


I 

r' 




{D^ + H^)dx + d^ + h‘^ 


-1 

Dl 


moE 


Ht\dx- 


mp 

Ae 


1 — 2Kp(j)p{l) 


M 


D^dx 


(3.19) 


2mo7eCi^ moKo(j)o{l){l + 2'^,fd‘^ ^ 


3d'^ 


-moKocj)o ( 1 ) 


(l + 7e)e ' (1 + 7 £)Mv^ 2(1 +7 s) ' M^/e' 

We set 1 — 2i^o0o(l)^ > | for 0 < £ < £ 2 , where £2 > 0 depending only on mp and (/>o(l)- As 

((1+27 max )^+ 3 )/i'o</'o(l) 


le 

W 


> 


4M 


for Theorem 11.11 the last three terms of (13.191) become nonpositive if 
holds true. Hence we obtain fll.32p by letting e' := min{£o+ 2 }- 

Now we claim that if Ip < for some 0 < 0 < 1, then d‘^{t) < I(t) < /(O) = Ip for all 
t > 0, where I{t) ;= | + H'^)dx + Notice that by fll.32p . 


A 

dt 


(D^ + H‘^)dx + d^ + h^ 




< 0 , 


(3.20) 


provided that d? < Ip < Assume Ip < for some 0 < 6* < 1. Then fl3.20p implies that 
d^{t) < lit) < /(O) = Ip < for all t > 0 (see Appendix III for the detail). Therefore, we 
conclude that | > |(1 — d) > 0 for all t > 0, and complete the proof of Theorem 11.31 

Proof of Corollary 11.41 
By Theorem 11.31 we have 


A 

dt 


+ H^)dx + d^ + h‘^ 


'-1 


< - 




D^dx 


1-1 


for t > 0 and 0 < E < e'. As for the proof of Corollary 11.21 we have 


d 

dt 


(D^ + H^)dx + d^ + h^ 




< -o' 


(D^ + H^)dx + d^ + h^ 


'-1 


for some o' > 0 depending only on 9 and the constants from the Poincare’s inequality. There¬ 
fore, we complete the proof of Corollary 11.41 


Appendix I 


Here we want to prove that ||hl||j ;^2 + is equivalent to ||h||^-i + ||p||j|^-i, which means 

that 

Cl (lIBII,,, + + ||fi||„-. < Cl (||B||,,1 + IliTUi) 

for some constants Cj > 0, j = 1,2 independent of D, H, n and p. Because h = 1(6 + fj) and 
p = |((5 — ?)), it is sufficient to show that 


Ci\\D\\^, < ||(5||^-i < ^ 2111111^2 and Ci||dd| 


Li ^ 


HZ 


-1 < CoWHl 


1,2 5 


(3.21) 
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for some constants Cj > 0, j = 1,2 independent of D, H, S and fj. For any fixed t > 0, we set 
S G H~^ ((—1,1)) and denote {S,v) and ||(5||j:^-i as follows: 


{S,v) = J Svdx for all v G {{—1,1)) , 


and 


" ll^ll i=ii_i 

where ( , ) denotes the pairing between and By the dehnition of D (see (ll.22p ). 
we have = 5, D{±l,t) = 0 and hence {S,v) = — Dvxdx for v G 1,1)) using 

integration by part. 

Now we claim that ||5||j;^-i is equivalent to ||Z1 ||l2 , which means that Ci||Zl ||^2 < ^ 

C' 2 ||/l||j ;^2 for some constants Cj > 0, j = 1, 2 independent of D and 5. For any v G H^{{—1, 1)) 

with \\ v \] h ^ = 1, 

/•i 


{6,v) 


Dvrdx 


'-1 


^ II-D II L2 11^x111,2 < ||hl||L2 . 


2 ||2 ii2 

Here we have used Holder’s inequality and the fact that 1 = Unllj^i = 11^11^2 + 11^3,11^2- Conse¬ 
quently, ||5||j^-i < ||/ 1 ||l 2 . On the other hand, let 


V{x,t) := [ D{s,t)ds - 

J~i 2 


^1 /•X 


D{s, t)dsdx , for x G (—1,1), t > 0 . 


'-1 j-i 


Then = D, Vdx = 0, and hence we have 


D^dx = 


VSdx 




'-1 


< IlilU-.IIPIU; 

— ^Il^lli?yiII-^IU e by Holder’s and Poincare’s inequalities, 

which implies ||T*||i ,2 < (711(511^:^-1. Note that = V^D, = 5, D{±l,t) = 0, and here we 
have used integration by parts for the first equality. Therefore, ||(5||j:^-i is equivalent to ||Zl||i 2 . 


Similarly, we may get the equivalence between 

dMH). 


H. 


-1 and ||Tf||L 2 and complete the proof of 


Appendix II 


Here we state the proof of Proposition 13.21 Multiply equations of n and p of fll.7p by 
n_ := min{?7,, 0} and p_ := min{p, 0}, respectively, and integrate them over (—1,1). Then 


1 d 

2 dt 


n dx = 


12 - fK^dx^ 

2 dtJ_J- 


rl 




{{n-)l - n-{n-)x4>x) dx, 
+ P-{P-)x(i)x) dx. 
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We use the interpolation inequality ||n||i ,3 < C'||m||]4^||m||]:^i for u G 1,1) to get 


n-{n-)a,(t)xdx 


'-1 


< 


i 


n_{n_)x(j)x\dx 

< II {ri_)x\\Ll{{-i,i))\\n- ||l3((_i,i))||0x||l6((_i,i)) 

< Cl||(?7,_)3;||i2((_l^l))||n_||^2((_l_l))||’^-||j^l((_l_l))|10x||L6((_l^l)) 

^ II 11^/^ II 11^/^ 11^ II 

< + C'2|h-|li2((_i^l))||0x||L6((_i4)), 

where Ci, C 2 depend only on the domain (—1,1). Using Poisson’s equation of 0 in fll.Tp with 
Robin boundary condition of fll.Sp . we have | 103 ;||l 6 ((-i,i)) < 1 ^ 3(1 + \W ~ PlUi((-i,i)))) where C 3 
depends only on the domain (—1,1), e, 7 ^, and 0o(±l)- Therefore, 


j < C'sC's (1 + ||n-p||i2((_i_i)))^y n^_dx 


(3.22) 


= f{t) / riidx, 


'-1 


where f(t) := C 2 C 3 (l + ||n — p||l 2 ((-i,i)))"^ ^ -^^((0) ^)) T > 0 since n,p G L°°(0, T; T^(- 

solve fll.Vll - fll.Sp . By fl3.22l) and the fact that n_{x, 0) = 0 for x G (—1,1), we have n_ = 0, i.e., 
n > 0. Similarly, we get 


1 d 

2 dt 


p^_dx < f{t) / ptdx, 


(3.23) 


'-1 


'-1 


and p_(x, 0) = 0. Therefore, p_ = 0, i.e., p > 0 and we may complete the proof of Proposi¬ 
tion [321 


Appendix III 

Here we prove that I (t) < I (0) = Iq < for f > 0 if / (0) = /q < holds true for some 

0 < 6^ < 1. Using fl3.20p . it is equivalent to show that ^ = [t > 0 : I {s) < for 0 < s < t} = 
[0, cxo) if / (0) = Iq < holds true for some 0 < 6^ < 1. By the continuity of function J, there 
exists fi > 0 such that / (t) < for 0 < t < ti, which satishes the condition of (I3.20p . Please 
note that (t) < I (t) for f > 0. Consequently, fl3.20l) implies that I {ti) < I (t) < I (0) < 
for 0 < f < ti- That is, ti G O'. Moreover, we claim that is open in [0,cxo). Suppose 

0 < t 2 ^ Then I (t) < fQ^- f < which implies [ 0 ,^ 2 ] C O'. By the continuity of 

function J, there exists d > 0 such that / (t) < for t 2 < t < t 2 + S, and the condition of 

fl3.20l) holds true for t 2 < t < t 2 + 6 . Hence by fl3.20l) . I (t 2 + S) < I (t) < I (^ 2 ) < for 

^2 < t < ^2 + d, which implies that [^ 2,^2 + 5) C O' and O' is open in [0, 00 ). On the other 
hand, it is trivial that O' is closed in [0, cxo) because of the continuity of function I. Therefore, 
0= = {f > 0 : J (s) < for all 0 < s < f} = [0, cxo) and fl3.20p gives I {t) < I (0) = /q < 
for f > 0. 


1 , 1 )) 
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